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We investigate the statistical mechanics of long developable ribbons of finite width and very small 
thickness. The constraint of isometric deformations in these ribbon-like structures that follows from 
the geometric separation of scales introduces a coupling between bending and torsional degrees of 
freedom. Using analytical techniques and Monte Carlo simulations, we find that the tangent-tangent 
correlation functions always exhibits an oscillatory decay at any finite temperature implying the 
existence of an underlying helical structure even in absence of a preferential zero-temperature twist. 
In addition the persistence length is found to be over three times larger than that of a wormlike 
chain having the same bending rigidity. Our results are applicable to many ribbon- like objects in 
polymer physics and nanoscience that cannot be described by the classical worm-like chain model. 



A shell or plate is a body whose two dimensions are com- 
parable to each other and much larger than the third. A 
rod, on the other hand, is a body whose two dimensions 
are comparable to each other and much smaller than the 
third. In each of these objects, this separation of scales 
leads to constraints and an energetic ordering of the dif- 
ferent possible deformations such as bending, twisting, 
stretching and shearing, and there are well accepted con- 
tinuum and statistical theories of their behavior. There 
is, finally, an intermediate class of objects where all three 
dimensions are widely separated from each other, which 
we call strips or ribbons. 

An example of such a chimeric object is a rectangu- 
lar sheet of paper of large aspect ratio (see Fig. [I]); its 
thickness is much smaller than its width which is itself 
much smaller than its length. The nanoscale world of- 
fers a number of compelling examples of these ribbon-like 
structures such as DNA, the secondary structures of pro- 
teins such as a-helices and /3-sheets, Carbon (graphene) 
and Molybdenum ribbons etc. A classical theoretical ap- 
proach to the study of these filamentous objects is em- 
bodied in the wormlike chain (WLC) model [1 , in which 
the polymer is assumed to have an inextensible back- 
bone, but with a flexibility governed by the persistence 
length £ p which expresses the exponential rate of decay 
of tangent-tangent correlation function. Its relative sim- 
plicity and precision in describing the entropic elasticity 
of many semiflexible polymers has made it very attrac- 
tive as an object of theoretical study. However, the ac- 
curacy of the WLC model for finite width biopolymers is 
still debated and a number of alternatives have been pro- 
posed to account for the torsional rigidity, the excluded 
volume and other features associated with a finite width 
[2]-[4]. In this Letter we consider a model for the statisti- 
cal mechanics of ribbon-like objects based on a classical 
geometric formulation of developable surfaces. In the 
limit when the thickness of a ribbon is much less than 
its width which is itself much less than its length, we can 
define a centerline. Then isometric deformations of the 
ribbon (i.e. those that preserve local distances so that 



the metric tensor is always Euclidean, gij = 5{j) can be 
described in terms of the orientation of the generators 
of the developable surface relative to the centerline [5] 
and account for the finite width of the ribbon. Because 
of its tensorial nature, developability strongly constrains 
the allowable deformations of a ribbon in a fundamen- 
tally different way than other milder constraints such as 
the "railway-track" models of polymers [3]. The most 
striking consequence of this, illustrated in Fig. [I] is that 
a developable ribbon can bend without twisting but can- 
not twist without bending, which as we shall see has far 
reaching consequences. 

To characterize a ribbon, we first consider the ge- 
ometry of its centerline in terms of the Frenet frame 
(t, n, b)(s) of the tangent, normal and binormal vector 
at each point along the curve parametrized by arc-length 
coordinate s [6]. This orthogonal triad evolves in space 
according to the Frenet-Serret equations: 

t' = kyi . n' = — ftt + rb, b' = rn . 

where k,t are the curvature and torsion of the space 
curve and the prime denotes the derivatives with respect 
to s. In physical ribbons, there is another distinguished 
frame, the material frame which is attached to each cross- 
section; the spatial rate of rotation of this material frame 
is known as the twist r e and is not the same as the tor- 




FIG. 1: (Color online) A twisted paper ribbon. In presence 
of twist the ribbon necessarily develop a non-zero curvature, 
as required by the Sadowsky functional. The Frenet frame 
(t, n, b) and the material frame (di, d2, ds) are sketched in 
red and green respectively. For a developable ribbon, the two 
frames coincide, so that the twist of the material frame and 
the torsion of the center line are identical. 
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sion of the geometric center line. For example, an arbi- 
trary twisted ribbon can have a straight center line, while 
a helically wound telephone cord may only have torsion 
and no twist. However, when the ribbon is deformed iso- 
metrically, the direction perpendicular to the centerline, 
which points along the binormal, also points along a ma- 
terial axis that corresponds locally to the only direction 
in which bending is possible, so that Frenet frame and 
the material frame coincide and r = r e . Equivalently, 
the projection of the curvature of the centerline on the 
material frame always vanishes in the infinitely stiff di- 
rection. 

The elastic energy due to isometric deformations of a 
naturally flat sheet of thickness h is given by the func- 
tional: F = \D J dA(Ki + ^2) 2 7 where K\ and ^2 are 
the principal curvatures, D = h 3 E/[12(l — v 2 )] is the 
bending rigidity, E the Young modulus of the material 
and v its Poisson ratio. The assumption of developabil- 
ity allows one to reduce the two dimensional integral to a 
one- dimensional integral along the centerline of the strip 
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where k and r are respectively the curvature and the 
torsion of the centerline, w is the lateral width of the 
ribbon and L the length of the centerline. In the limit 
this yields the Sadowsky functional [5]: 
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As anticipated, the nonlinear coupling between curvature 
and torsion in implies that an inextensible ribbon can 
bend without twisting (r = 0, n 7^ 0), but cannot twist 
without bending (r 7^ —> K ^ 0). 

To make the Sadowsky's functional amenable to 
numerical simulation of a discretized analog, we con- 
sider a chain of N = L/a consecutive vertices x = 
{xi, X2 . . . xat} separated by a fixed length a. At each 
point x^ in the chain, the discrete Frenet frame is given 
by: = (t^, n^, hi), are orthogonal 3x3 matrices 
whose columns are given by the column vectors t^, and 
b^. Calling 0i and fa the angles formed by consecutive 
tangent and binormal vectors (i.e. t{ • t^_i = cos^ and 
b^ • b^_i = cos^), the Frenet equations for the moving 
trihedron are expressed as F^ = F^iR^ [9 , where 



cos 0j — sin 0j 
Hi = sin 0i cos fa cos 0i cos c 
sin 6i sin fa cos 0i sin ^ 




(2) 



We note that ti • n^_i = sin^cos^ and ti • b^_i = 
sin 0i sin fa. Thus 0i and fa are also the polar and az- 
imuthal angle of a system of spherical coordinates (see 
inset of Fig. [2j. The discrete analog of curvature and tor- 




FIG. 2: (Color online) The quantity A denned in the text as 
a function of /?, calculated from the exact integral expression 
given on the last page. With exception for the two limits f5 — > 
0, 00, A is always negative and the associated eigenvalues 
Ai,2 = \(b =b a/A) are complex and conjugate. In the inset 
the discrete Frenet frame. 

cos^) and rf = a _2 |bi - b^_i| 2 = 2a _2 (l - cos< 
energy can thus be written as: 
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where f3 = Dw / '(afc^T). This serves as a starting point 
for Monte Carlo simulations of a chain composed of N 
segments of unit length and its conformation is updated 
using pivot moves through a Metropolis algorithm. For 
every chain we perform 10 6 Monte Carlo sweeps, the first 
half of which is devoted to equilibration. Fig. |3]shows the 
tangent- tangent and binormal-binormal correlation for a 
chain of N = 100 segments at various temperatures. We 
see that while (b n -bo) always exhibits exponential decay, 
(t n • t ) has oscillations superimposed on an exponential 
decay. Upon raising the temperature, the frequency in- 
creases, but the exponential damping becomes stronger. 

To understand this result, we first derive the 
tangent- tangent and binormal-binormal correlation func- 
tions from the iterated transfer matrix: (t n - to) = 
(R!R 2 . . . Rn)n and <b n • b ) = (R1R2 . . . Rn)s3 H, 
noting that since the energy ([3| is a sum over individ- 
ual segments, the Boltzmann measure exp(— F/ksT) is 
factorizable and (R1R2 . . . R n ) = (R) n . Moreover, for 
sufficiently large spatial separation the behavior of the 
correlation functions will be dictated by the largest eigen- 
value in the spectrum of R. Assuming (sin fa) = 0, the 
transfer matrix R becomes block- diagonal and the eigen- 
values are given by the roots of the characteristic equa- 
tion A 2 — bX + c = with b = (cos 6) + (cos 6 cos fa) and 
c = (cos 0) (cos cos fa) + (sin#)(sin#cos(/>), so that 
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where A = b 2 — 4c. Thus the binormal-binormal correla- 
tion function always exhibits exponential decay: 
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FIG. 3: (Color online) Tangent-tangent (left) and binormal- 
binormal (right) correlation functions from Monte Carlo sim- 
ulations of a polymer chain of N — 100 segments at T — 0.1 
(black), 0.2 (blue), 0.3 (red) and 0.4 (green). The graphs of 
(t n • t ) have been translated along the vertical axis to avoid 
overlaps. 



where £~ l = —a -1 log (cos <j>) is the torsional persistence 
length and s = na. On the other hand, the sign of A 
dictates the behavior of the tangent-tangent correlation 
function. When A > 0, the eigenvalues are real and 
(t n • t ) has the typical exponential decay of a WLC. 
When A < the eigenvalues are complex and (t n • to) 
has the form: 
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with k = a 1 atan(V— A/6) and £ p 
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exact numerical calculation of A, shown in Fig. [2] as 
function of /?, shows it is always negative so that the 
tangent-tangent correlation function is always oscillatory 
at any non-zero temperature. 

This oscillatory behavior is reminiscent of spin chains 
with competing ferromagnetic and antiferromagnetic in- 
teractions. In these systems the oscillatory behavior of 
the two-point correlation functions denotes the existence 
of a helical phase in which the magnetization varies pe- 
riodically in space. The critical point that divide the 
disorder phase and the helical phase is traditionally re- 
ferred to as "Lifshitz point" [ID] . Analogously, an oscil- 
lating tangent-tangent correlation function in our poly- 
mer model implies the presence of an underlying helical 
structure that persists at any finite temperature even in 
absence of a spontaneous twist, i.e. so that the Lifshitz 
point here is T = oo. In Fig. [4] we show that our numer- 
ical simulations capture this feature, due to the strong 
geometrical constraint of develop ability, i.e. the ribbon 
is isometric to a flat strip at any temperature, and thus 
bends when twisted. The notion of a Lifshitz point has 
been used earlier to described the low temperature phase 
of a "railway-track" polymer consisting of two parallel 
WLCs connected by a rigid bond [3], and is marked by an 
oscillatory binormal-binormal correlation function which 
results in a "rod-kink" structure in which the ribbon is 
at every point along its contour either twisted and un- 
bent (rod) or bent and untwisted (kink). This is quali- 
tatively different from our results where it is impossible 
for the ribbon to be twisted and unbent simultaneously, 
due to the fact that the constraint of developability is 






FIG. 4: (Color online) The typical conformation of a ribbon- 
like polymer after 10 6 Monte Carlo iterations. (Left) = 
0.1 and (right) = 1. At low temperature the ribbon 

fluctuates around an underlying helical shape. The radius 
and the pitch of the helix are set by the persistence length £ p 
and the wave- number k and are proportional to (3a and f3 1 ^ 2 a 
respectively. 



implemented exactly in Eq. ([I]), unlike in any previous 
attempts. 

Next we give an analytical estimate of the three rel- 
evant quantities £ p , £ T and k in the low temperature 
regime. The partition function for the ribbon is given 
by Z = with: 
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with d£i = d(j)d0 sin the infinitesimal solid angle. The 
analytical calculation of thermal average is complicated 
by the coupling term between curvature and torsion, but 
this can be removed by rewriting the r 4 /^ 2 term as a 
Gaussian integral of the form: 
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With this substitution, curvature and torsion are now 
decoupled and the single-node partition function reads: 
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with A = (1 + y 2 ) and \i = 2(1 + iy), so that the rib- 
bon is equivalent to a wormlike polymer with finite but 
complex- valued bending and torsional stiffness. Calcu- 
lating the solid-angle integral in the previous expression 
yields 
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where Io and 7 are respectively the modified Bessel func- 
tion of zeroth order and the incomplete gamma function 
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[TT] . Analogously one can calculate the ensemble average 
( ' ) = z i X I If ( • ) e-^ a2 ^ 2+2r2+r4 /* 2 ) which yields: 
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(sin0cosm</>) = (^) V* dye-^I m (^) / 2 (/?A) , 
where m is an integer in {0, 1} and 
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with erfi the imaginary error function [11 . At low tem- 
peratures (/? ^> 1), asymptotic expansions of the func- 
tions Iq, 7 and erfi yield 



(cos 6 cos mq 



(sin cos mq 



21(1 -4m 2 ) 4(7- 173m 2 ) 



16/3 

128/3[40/3 - 7(2 ■ 



5 + 128/3 ' 
-99(1 -4m 2 ) 



12v^/3 3 / 2 (5 + 128/3) 
from which we find the persistence lengths £ p and t T 
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A plot of £ p , £ T and /c versus temperature is shown in Fig. 
[5] We see that as a consequence of the geometric coupling 
between curvature and torsion in a developable ribbon, 
the persistence length is a factor (|§ — §f ) = 3.476 
larger than that of a WLC at the same temperature. 

Our study has a set of clear experimental predictions 
for the equilibrium statistical mechanics of ribbon-like 
polymers where all three length scales are well separated. 
Firstly, the tangent- tangent correlation function exhibits 
an oscillatory decay at any finite temperature, implying 
an underlying helical structure even in absence of any 
natural twist at zero temperature. For a ribbon with fi- 
nite extensibility, even if the helical phase disappears at 
high enough temperatures, the Lifshitz point should per- 
sist but is shifted to T = T* , and the mean field estimate 
for the oscillatory wave number k ~ [kB(T—T*)/Dwa] 1 / 2 
The helical structures that arise herein are reminiscent 
of those found in biopolymers and proteins, and suggests 
a natural generalization of our theory to include short 
range forces that will stabilize these structures. In the 
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FIG. 5: (Color online) Persistence length and torsional persis- 
tence length (left) as a function of (in units of a). (Right) 
Tangent- tangent correlation function wave-number k (in units 
of a -1 ) as a function of Dots corresponds to numerical 

data and lines to analytical predictions for the low temper- 
ature regime. For biopolymers of Young modulus E « 100 
MPa, h ~ 1 nm and w/a ~ 1, the value of /3 _1 in the plots 
spans the temperature range T G [0, 10 3 ] K. 



presence of forces, the geometrical structure of ribbons is 
expected to lead to qualitatively differences in the force- 
extension curves of these ribbons as well as in their mo- 
tion, topics of future study. 
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